PROCEEDINGS 

0F SCIENCE 



Small, dark, and heavy: But is it a black hole? 



Matt Visser* 1 

School of Mathematics, Statistics, and Operations Research, Victoria University of Wellington, 
New Zealand 
E-mail: matt . visser@msor . vuw . ac . nz 

o 

CN Carlos Barcelo 

C*) Instituto de Astrofisica de Andalucia, CSIC, Camino Bajo de Huetor 50, 18008 Granada, Spain 

1^ £-mai/: carlos@iaa.es 



Stefano Liberati 

International School for Advanced Studies (SISSA), Via Beirut 2-4, 34014 Trieste, Italy 
and INFN, Sezione di Trieste 



1 February 2009; 10 February 2009; WT^C-ed February 9, 2009 



E-mail: liberati@sissa . it 

6JQ Sebastiano Sonego 

Universita di Udine, Via delle Scienze 208, 33100 Udine, Italy 
£N) E-mail: sebastiano.sonego@uniud.it 

> 

Astronomers have certainly observed things that are small, dark, and heavy. But are these objects 
really black holes in the sense of general relativity? The consensus opinion is simply "yes", and 
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1. Introduction 

Do alternatives to standard classical black holes exist? Can one "mimic" a black hole with arbitrary 
accuracy? There is a rather limited set of (arguably) viable alternatives: 

• Quark stars [1], Q-balls [2], strange stars [3]. 

• Boson-stars [4]. 

• Gravastars: Mazur-Mottola variant [5], and Laughlin et al. variant [6]. 

• Fuzz-balls: Mathur et al. variant [7], and Amati variant [8]. (See [9] for a survey.) 

• Dark stars/quasi-black holes [10, 11]. (For related ideas, see [12, 13]). 

While there are close inter-relationships between these various models, in this article we will specif- 
ically focus on our own proposal [10, 11], and give an informal overview of the situation. We shall 
re-assess and (hopefully) re-invigorate an old line of argument: What effect does quantum physics 
have on the collapse of a classical star? Is semi-classical collapse [10, 12, 13, 14, 15] qualitatively 
different from classical collapse [16]? 
In general we can certainly write 



and separate the expectation value of the stress-energy-momentum tensor into a contribution from 
some suitably chosen vacuum state, plus a contribution from the excitations above that vacuum 
state. For instance, for an uncollapsed star 



The vacuum polarization effect (0| t a b | 0) is utterly negligible in an ordinary uncollapsed star. (This, 
after all, is why we can get away with just solving the classical Einstein equations most of the 
time.) Does this remain true during collapse? Even if the vacuum polarization does remain small, 
it might still have a significant effect on the location and/or existence of event horizons [14, 15]. 

Now this point of view, while certainly historically respectable, does deviate significantly 
from the present "consensus opinion", at least in the general relativity community, so before one 
gets started there are a number of preliminary issues that should be dealt with. 

2. The Fulling-Sweeny-Wald (no-singularity) theorem 

There is a widespread feeling in the general relativity community that semiclassical quantum 
back-reaction effects are always small, and never enough to significantly alter the classical picture 
of collapse to a black hole. (See figure 1 for an appropriate Carter-Penrose diagram.) When 
pushed, members of this community ultimately point to the Fulling-Sweeny-Wald no-singularity 
theorem [17] as the basis for this assertion. 

Phrasing the Fulling-Sweeny-Wald theorem rather loosely: "In quantum field theory on a 
curved spacetime everything is finite at an event horizon, and all the way down to either a singu- 
larity or a Cauchy horizon." The technical content of this theorem is based on showing that the 



G ab = S7t(Y\f ab | y) 



(1.1) 




(1.2) 
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Singularity 




Figure 1: Standard Carter-Penrose diagram for an astrophysical black hole formed via stellar collapse. 

Hadamard form of the two-point function for a quantum field is not affected by the presence of an 
event horizon. (So for a Hadamard quantum state, the renormalized stress-energy tensor [RSET] is 
automatically finite.) Unfortunately, for the question we want to raise, the Fulling-Sweeny-Wald 
theorem also "begs the question", for what it shows is that if an event horizon forms, then the quan- 
tum field theory is well behaved there. But this is not the same as showing that an event horizon 
will naturally form in semiclassical collapse: A finite but very large contribution from the vac- 
uum polarization term (0|r a ^|0) in equation (1.2), while perfectly in agreement with the theorem, 
would significantly alter the dynamics of collapse. In fact, compact horizonless objects, and/or 
naked singularities, are also fully compatible with both the hypotheses and the conclusions of the 
Fulling-Sweeny-Wald theorem. 

3. Horizons and the choice of the quantum vacuum state 

Our specific proposal for a black hole mimic is ultimately related to the appropriate choice for 
the quantum vacuum state |0). Common candidates states are: 

Boulware vacuum: This is singular at any Killing horizon. The renormalized stress-energy- 
momentum diverges as r — > 2m. 

Unruh vacuum: This is designed to be well-behaved at any future Killing horizon, so in particular 
the renormalized stress-energy-momentum is finite there. 

Neither of these quantum states contain particles in the vicinity of past null infinity . Accord- 
ingly, the key constraint that we shall adopt for our global vacuum state is that it contain no particles 
in the vicinity of . Additionally, far in the past, near past timelike infinity i~, when spacetime is 
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static (and nearly Minkowskian), the quantum state should exhibit properties qualitatively similar to 
those of the Boulware vacuum. This demand is physically appropriate, since the Boulware vacuum 
is the one that best describes physics in the presence of a static self-gravitating object. Note that if 
we were instead to choose a state that in the asymptotic future behaves like the Unruh vacuum, then 
this would presuppose the formation of a horizon, which however is precisely the issue we wish to 
investigate. Such a choice would anyway amount to making a teleological statement. Apart from 
these particular issues, there is an increasing consensus, or at the very least a suspicion, within 
the general relativity community that event horizons are simply the wrong thing to be looking 
at. Sometimes apparent horizons [18, 19, 20] (or better yet, dynamical horizons [21], or trapping 
horizons [22]) are better candidates for characterizing the black hole. (See also [23, 24, 25].) 



4. Our specific calculation 

The metric for the spacetime of a spherically symmetric collapsing body can be written in 
Schwarzschild coordinates as 

g=- e - 2 ^ (l-2m(r,t)/r)dt 2 + ^ ^ . . +r 2 (d0 2 + sin 2 0d(p 2 ) . (4.1) 

\—2m(r,t)/r 

In these coordinates an apparent/dynamical/trapping horizon, if it forms, is characterized by 

2m(r,t) 

— ^ = 1 . (4.2) 
r 

In contrast, an event horizon, if it forms, can only be found by back-tracking from future null 
infinity, J? + . 

In the standard classical conformal Carter-Penrose diagram for the collapse process, figure 2, 
one truncates the diagram at the centre, r = 0, and modes of any field residing on the spacetime 
are said to "bounce" off the centre. A central part of the analysis is then to relate the affine null 
parameter W on to the affine null parameter u on via some function W = p(u), which thus 
encodes a good fraction of all the physics. 

For technical reasons we prefer to work with a more "symmetric" version of the Carter- 
Penrose causal diagram. In this version, figure 3, modes propagate straight through the centre 
of the collapsing star (located at r = 0), so that J^~ ft is connected to =^t ht , and vice versa. Affine 
null coordinates on are now related to affine null coordinates on by: 

U = p{u) ; W = p{w) . (4.3) 

Furthermore, we shall work in a 1 + 1 dimensional truncation of 3 + 1 dimensional spacetime in 
order to keep the calculation tractable. It will then be convenient to extend the coordinate r in such 
a way that on the right-hand part of the diagram r > 0, while in the left-hand part r < 0. 

For a massless real scalar field, one can expand the field operator using a set of modes that 
near J^ left take the form 
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Figure 2: Standard Carter-Penrose diagram for the collapse process: Modes "bounce" off the centre of the 
star. We deliberately leave the "late time" portion of the diagram vague and ambiguous — since we do not 
want to pre-judge what the system settles down to. 



with Q. > 0; these are appropriate to define particles in the asymptotic past, before collapse takes 
place. Neglecting backscattering, such modes take, near J'^ghv trie f° rm 

^'^ wmm^r ^- (4 - 5) 

However, near J^ ht , the modes appropriate for defining particles are not the cp^, but others that 
we denote Ya, with the asymptotic form 

""»<'•" " wr-Wnr """" ■ <4 - 6) 

Hence, the state which does not contain particles on (defined using the modes (pa) turns out to 
contain particles on J? + (defined using the Y®)i provided that p(u) is a non-trivial function, such 
that a (pa mode contains negative-frequency contributions when Fourier-analysed in terms of the 
Ya ■ Defining the w-dependent frequency on , 

<o(u,Q) = p(u) O, (4.7) 

associated with a q>Q mode, one can see that mode excitation takes place provided that the adiabatic 
condition, 

1 v ; 71 < i, (4.8) 
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Figure 3: Our preferred Carter-Penrose diagram for the collapse process: Modes propagate straight through 
the centre of the star. We deliberately leave the "late time" portion of the diagram vague and ambiguous — 
since we do not want to pre-judge what the system settles down to. U-W coordinates are best in the bottom 
corner of the diagram, in fact near all of . In contrast, u-W coordinates are best in the upper-right region, 
near (And U-w coordinates are best in the upper-left region, near J^ja.) 



is violated. This occurs for frequencies smaller than 

BoM-M. (4.9) 
PW 

which can then be thought of as a frequency marking, at each instant of retarded time u, the sep- 
aration between the modes that have been excited (those with £2 <C £2o) an d those that are still 
unexcited (£2 3> Oo). Intuitively, one may think that there is still an infinite "reservoir" of high- 
energy Boulware-like modes, and that if mode excitation is not sufficiently rapid they will make 
a potential obstruction to horizon formation. Indeed, calculations in static models [26] show that 
such modes lead to an energy condition-violating renormalised stress-energy-momentum tensor 
that diverges as one approaches the horizon. We feel however, that this is far too naive a picture, 
being based on results obtained in static spacetimes, and that it can be taken at best as a hint that 
one should carefully check how the renormalised stress-energy-momentum tensor behaves when 
the horizon is just about to form. So, let us now turn to a calculation that takes dynamics explicitly 
into account. 

The spacetime metric can be written using either the set of coordinates (U,W), or (u,W): 

g= -C(U,W)dUdW = -C(u,W)dudW . (4.10) 
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This gives 

C(U,W) = ^1, (4.11) 
p{u) 

where, for events lying outside the collapsing star, C(u,W) is the metric coefficient of a static 
spacetime. For any massless quantum field, the renormalised stress-energy-momentum tensor cor- 
responding to a quantum state that behaves like the Boulware vacuum asymptotically in the past 
has components [26] 

T uv -C 1 / 2 dlC- 1 ' 2 , T ww ~C [ l 2 dlC- l l 2 , T uw ocCR, (4.12) 

where R is the curvature scalar, and the numerical coefficients depend on the specific type of field 
being considered. 

The component with the most interesting structure is Tuu- Using equation (4.11) and the 
property du = p~ Y d u one finds 

C x l 2 d 2 C- x ' 2 = j 2 (c l l 2 d 2 C- l l 2 -p 1 ' 2 d 2 u p- 1 ' 2 ) . (4.13) 

The key point is that the first term within brackets on the right hand side of equation (4. 1 3) is a static 
contribution due to the Boulware-like modes, while the second one arises because of the dynamics 
of collapse. These two terms, separately, would lead to an arbitrarily large Tuu as the horizon is 
approached, because p tends to vanish in this limit. However, if (and only if) the horizon forms, 
then the leading contributions of C 1 ' 2 d 2 C~ 1 ' 2 and p l ^ 2 d 2 p~ l l 2 exactly cancel against each other. 

For the computation, it is convenient to work in a chart that is regular at the horizon (if it 
forms), so that the regularity of the stress-energy-momentum tensor can be inferred just by the 
finiteness of its components. This is the case for the Painleve-Gullstrand coordinates (x,f), in 



terms of which the metric is [27, 28, 29]: 

g= -c(x,t) 2 dt 2 + (<hc-v(x,t)dt) 2 . (4.14) 

A rather technical computation gives [10]: 

T, t = U 2 Tuu + 2U t W t T uw + W 2 Tww (4.15) 

= {c + v) 2 U 2 Tuu-2{c 2 -v 2 )U x W x Tuw + {c-v) 2 W 2 Tww (4.16) 

= p 2 Tuu-2pTuw + T W w ; (4.17) 

T tx = U, U x Tuu + (U t W x + U x W t ) T uw + W, W x T ww (4. 1 8) 

= -(c + v)U 2 Tuu-2vU x W x T uw + (c-v)W 2 Tww (4.19) 

p 2 2pv 1 

= ; — Tuu + -j j Tuw H Tww ; (4.20) 

c + v c z — v z c — v 

T xx = U 2 Tuu + 2 U x W x T uw + W 2 Tww (4.21) 

p 2 p 1 

= -; — i — ry T V u - 2 -s j T V w + 7 T W w ■ (4.22) 

(c + v) z c z — v z (c — v) z 
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Since at a horizon c + v — > and c — v — > 2c, the term that presents, potentially, the highest degree 
of divergence is the one proportional to Tyjj in T xx . The other potentially dangerous coefficients 
p~ 2 and p~ x that appear in Tuu an d Tyw are cancelled by corresponding factors in the expres- 
sions (4.17), (4.20), and (4.22). 

In the rest of the calculations, we assume that c(x) = 1 and place the horizon (when it exists) 
at x = 0. Then, assuming that a horizon indeed forms, we can expand v(x) as 



V (jc) = -l + KX+ff(x 2 ) , 



(4.23) 



where K can be identified with the surface gravity [28, 29]. The static contribution in equation 4.13 
is 



(4.24) 



and this, taken alone, would cause the T tx and T xx coefficients to diverge. However, under the 
hypothesis of horizon formation one also has [10] 



p(u) = U B -Ai &- Ku + ^&- lKU + A3 ~- 3Kl < 



+ > 



2 3! 

where Uu, A\ > 0, Ai, A3 are constants. The dynamical term in equation (4. 13) is then 



(4.25) 



+ 



IA3 3 /A2V 
2Ai 4 UJ 



K 2 e -2m + 0, e -3Ku 



(4.26) 



Replacing the expressions (4.24) and (4.26) into equation (4.13), one sees that the dominant terms 
K 2 /4 cancel against each other, and one remains with a finite contribution that depends on the 
details of collapse: 



so that 



passuming horizon formation 



■^assuming horizon formation ( 



1 



1 A 3 3 
"2 A7 + 4 V Ai 



1A 3 3 (A 2 
2 Ai 4 1 Ai 



K 2 e- 2Kl ' + . 



e -2Ku + 



(4.27) 



(4.28) 



with these expressions holding outside the surface of the collapsing star. Furthermore, it is easy to 
realise that this contribution is inversely proportional to the square of the speed at which the col- 
lapsing body crosses its gravitational radius. Hence, for a very slow collapse there is a realistic and 
concrete possibility that the (energy-condition-violating) renormalised stress-energy-momentum 
tensor, although finite, could lead to significant deviations from classical collapse when a trapping 
horizon is just about to form. 

In order to reinforce this claim, let us consider a case in which the horizon never forms at any 
finite time, but is only approached asymptotically in the limit t — > +00. In particular, we shall be 
interested in an exponential approach [29], where the radius of the star depends on time as 



r(t) = 2M + B^ KDt 



(4.29) 
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with B and Kb positive constants. After a brief calculation [10, 29] this leads to 

p(u) = U u -A ie -^ u , (4.30) 

where 

K en =^<K. (4.31) 

Here K^ff can be thought of as a "reduced surface gravity". Interestingly, although no true horizon 
ever forms, one still gets a Hawking-like flux of Planckian radiation at the temperature [29] 

Keff KKn 

T = — = - (4 32) 

271 27t(k-+k: d ) 

Of course, outside the star, the calculation of C 1 / 2 d 2 Cr x l 2 again yields the same result as in 
equation (4.24). However, for the second contribution we now have 

K 2 

p l/2 dl p- x ' 2 = ^ + 0(e- 2K ^«) , (4.33) 

so there is no longer a perfect cancellation between the dominant terms in the static and dynam- 
ical contributions. Indeed, at the leading order, the renormalised stress-energy-momentum tensor 
outside the star in the limit x — > (that is, r — > 2M, or t — > +oo) behaves as 



wu 



K 2 -K 2 



4p 2 

so that 



(4.34) 



4k 2 x 2 

Note that this result is not in contradiction with the Fulling-Sweeny-Wald theorem [17], because a 
strict divergence appears only for t = +oo, i.e., at the boundary of spacetime. However, the renor- 
malised stress-energy-momentum tensor gives an arbitrarily large (albeit finite) energy-condition- 
violating contribution to the right hand side of the semiclassical Einstein equations as the horizon 
formation condition 2M/r = 1 is approached. 



5. Implications and Discussion 

In the standard collapse, you can argue that the renormalised stress-energy-momentum tensor at 
horizon-crossing, as felt by infalling matter, is negligible provided: 

1. the quantum state is of the Hadamard form (which we have by assertion); 

2. matter is basically freely-falling; 

3. the equivalence principle holds. 

The first point tells you that the quantum vacuum has the same ultraviolet form as in Minkowski 
spacetime, the second point tells you that matter is approximately in a locally inertial frame, and 
the third point tells you that the local renormalised stress-energy-momentum tensor the matter then 
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Figure 4: Standard Carter-Penrose diagram for an evaporating black hole. Don't bother asking what hap- 
pens at the endpoint of the evaporation process — in the standard causal picture there is no definite answer. 

"feels" must be approximately the same as in Minkowski spacetime — i.e., approximately zero 
(after renormalization). 

In contrast, our result is saying that large deviations from this standard conclusion can arise 
if matter is not freely falling, but is significantly accelerated (as, by self-consistency, it must be 
to sustain itself against the gravitational attraction). So we are explicitly violating point 2, (while 
we explicitly keep point 1, and implicitly keep point 3). So if the surface of the star deviates 
significantly from free-fall, then a large stress-energy-momentum builds up, which can force it 
further away from free-fall — either stopping or exponentially delaying the collapse. Precisely 
predicting what happens in a specific collapse scenario relies on extremely messy model-dependent 
physics. 

Indeed our calculation seems to suggest that if, during the late stages of the collapse, matter is 
far from free-fall, then a growing RSET can lead to a late time collapse history very different from 
the classically expected one, possibly leading to a form of asymptotic collapse of the type suggested 
in (4.29). It might even be that this is the solution preferred by nature; this might be due to new 
particle physics effect coming into place in the late stages of most stellar collapses. In this case the 
conformal diagram describing the gravitational collapse scenario would not be the standard one of 
figure 4, but rather that reported in figure 5. This object would then be a "quasi-black hole" (not 
to be confused with the homonymous objects proposed by Lemos and Zaslavskii [30], which are 
static solutions), an object which would not only closely mimic the classical geometry of a black 
hole, but also, (if the collapse law is exponential at late times), mimic its quantum effects such as 
Hawking radiation. 

To place our results in a broader perspective: Many physicists are now (for numerous indepen- 
dent reasons) arguing against the standard Carter-Penrose diagram, figure 4, for the formation and 
evaporation of a semi-classical black hole. Apart from the "physics challenged", (whom we shall 
quietly discount), there are hints from analogue spacetimes [29], from loop quantum gravity [31], 
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Figure 5: Conformal diagram of the spacetime for a quasi-black hole. The solid line represents the sur- 
face of the collapsing object; the dotted line is at r = 2M(t), where M(t) is the instantaneous mass of the 
object as measured from J !+ \ dashed lines correspond to (Schwarzschild) t = const hypersurfaces. The 
period of evaporation appears short because of a distortion induced by the representation, but actually cor- 
responds to a very long lapse of time, as one can see from the fact that the lines at t = const crowd around 
it. This diagram, while nonstandard, is nevertheless compatible with current astrophysical observations of 
gravitationally active collapse products. 

from string-inspired models [32], from attempts at unitarity preservation in our own domain of 
outer communication [33], from one-loop curved-space quantum field theory [34], and from ab- 
stract studies of the nature of horizons [35], all hinting at a more subtle history for collapse and 
evolution. (Canonical versions of alternative causal structures are given by the Carter-Penrose 
diagrams of figures 6 and 7, and the double-null diagram of figure 8.) Unfortunately, when at- 
tempting to move beyond qualitative statements of this type, specific predictions are frustratingly 
model-dependent, but there is some "wriggle room" for interesting new physics. 

On a cautionary note, we should point out that several authors have looked at the question of 
what observational signals for black hole mimics might look like [36, 37]. Critically, once you 
add rotation, the ergosurface is probably more important than the "would-be horizon". There is the 
very real risk of significant ergoregion instabilities [38]. 

In summary, what our calculation suggests is that it might be possible to have a black hole 
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Figure 6: Ashtekar-Bojowald version of the Carter-Penrose diagram for an evaporating black hole. The 
shaded region represents a region of Planck-scale curvature, and possibly large metric fluctuations. 



without having a black hole — a configuration that is a black hole for (almost) all practical pur- 
poses, but might be missing the one key ingredient of having a horizon. Deep issues of principle 
remain, and it will be very interesting to see how the whole area of black hole mimics develops 
over the next few years. 
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